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We explore the pairing properties of the one-dimensional attractive Hubbard model in the pres- 
ence of finite spin polarization. The correlation exponents for the most important fluctuations 
are determined as a function of the density and the polarization. We find that in a system with 
spin population imbalance, Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-type pairing at wavevector 
Q — |fcp,t — foul is always dominant and there is no Chandrasekhar-Clogston limit. We then in- 
vestigate the case of weakly coupled ID systems and determine the region of stability of the ID 
FFLO phase. This picture is corroborated by density-matrix-renormalization-group (DMRG) sim- 
ulations of the spatial noise correlations in uniform and trapped systems, unambiguously revealing 
the presence of fermion pairs with nonzero momentum Q. This opens up an interesting possibility 
for experimental studies of FFLO states. 
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I. INTRODUCTION 

In two-component Fermi systems with attractive in- 
teractions, in which both species share the same Fermi 
surface, pair formation is successfully described by the 
Bardeen- Cooper- SchriefTer (BCS) theory [1]. This sce- 
nario is, for instance, realized in conventional supercon- 
ductors, where electrons of opposite spin and momentum 
form Cooper pairs. In contrast, if the Fermi surfaces 
are different, the situation is subtler because not all of 
the fermions can form pairs. The question thus arises 
as to whether a superfluid phase exists in this setup 
and if so, what its properties are. On the theoretical 
side, different answers to this question have been pro- 
posed [2, 3, 4, 5, 6, 7, 8] ranging from exotic superfluid 
phases to more conventional phase separations. Partic- 
ular attention has been devoted to the so-called Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO, sometimes also de- 
noted LOFF in the literature) state, in which fermion 
pairs with nonzero momentum form an inhomogeneous 
superfluid phase. Initially developed for superconductors 
in magnetic fields, this theory has been applied to heavy 
fermion systems [9, 10] and dense quark matter [11]. 

Recent advances in methods for trapping and control- 
ling ultracold atoms have opened up the possibility of ex- 
perimentally observing population imbalanced mixtures 
of ultracold fermions [12, 13, 14, 15, 16]. Compared to 
solid state materials, cold gases allow an unprecedented 
control over the Fermi surface mismatch by adjusting the 
population imbalance 

_ Nj — Ni 
P iVf + iVj 

Here N a denotes the number of particles in one of two hy- 
perflne states labeled by a pseudo-spin index a G {|, |}. 
So far, these experiments have not revealed any of the 
proposed exotic superfluids, but have concluded that the 



system either phase separates, forming an unpaired nor- 
mal state immersed in a BCS superfluid [13], or forms 
pairs that do not condense [16]. Although theoretical 
studies of the stability and extent of the FFLO state 
in the phase diagram have not yet reached a consen- 
sus, it is clear by now that it is best realized at low 
temperatures and in quasi one-dimensional (ID) geome- 
tries [17, 18, 19, 20]. Given the disagreement between 
theoretical predictions and the difficulty of unambigu- 
ously identifying exotic phases encountered in experi- 
ments, we believe it lends insight to study the simplest 
possible system that exhibits FFLO-type pairing: the 
ID Hubbard model with attractive interactions. On 
the one hand, this system is well suited for theoreti- 
cal studies because it is exactly solvable via the Bethe- 
Ansatz and is numerically exactly treatable in density- 
matrix-renormalization-group (DMRG) simulations. On 
the other hand, it has now become possible to realize it 
experimentally [21]. This paper consists of two parts: 
In the first, we derive various correlation exponents as a 
function of the density and the polarization and demon- 
strate that the superfluid FFLO correlations dominate 
at all polarizations and that there is no Chandrasekhar- 
Clogston limit [22, 23] in the attractive ID Hubbard 
model. We also address the stability of the FFLO phase 
in weakly coupled chains. In the second part, we cal- 
culate the experimentally accessible spatial noise corre- 
lations quantitatively for uniform and trapped systems 
and show that FFLO states leave a distinct fingerprint, 
making it an ideal tool to detect and characterize exotic 
pairing phases [24, 25]. 



II. THE MODEL 

We concern ourselves with a population-imbalanced 
ID Fermi gas with attractive interactions, which can be 
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TABLE I: Correspondence between the correlation functions of the attractive and the repulsive Hubbard model. Note that the 
mapping holds for all three types of the triplet superfluidity. From Ref. 26. 



Attractive interaction 


U < 


U>0 


Repulsive interaction 


density correlations (CDW) 
singlet superfluidity (SS) 
triplet superfluidity (TS) 
longitudinal spin correlations (1SDW) 
transverse spin correlations (tSDW) 
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longitudinal spin correlations (1SDW) 
transverse spin correlations (tSDW) 
triplet superfluidity (TS) 
density correlations (CDW) 
singlet superfluidity (SS) 



described by the Hubbard model 

H = -* E { c L c i+i,« + Hc ) + u E n <>t n u • (!) 

1,(7 i 

Here cj a (c t a ) creates (destroys) a fermion with pseudo- 
spin a on site / and ni y(T = c\ a c { a is the occupation num- 
ber operator. The parameters t and U < characterize 
the nearest-neighbor hopping and the attractive on-site 
interaction, respectively. Furthermore, we specify the po- 
larization p and the particle density n = (iVj + N\)/L, 
L being the length of the chain. The physics of the 
attractive Hubbard chain for zero polarization is quite 
transparent: the on-site attraction U favors locally dou- 
bly occupied sites. These doubly occupied sites can lead 
to 2k f charge-density wave (CDW) fluctuations and to 
singlet superfluid (SS) fluctuations of pairs with zero mo- 
mentum. Away from half filling, the pairing fluctuations 
dominate for all densities when U < 0. Since spin excita- 
tions are gapped in this phase, spin density wave (SDW) 
and triplet superfluid (TS) correlations decay exponen- 
tially for U < 0. The correlation functions are defined 
in Table I. As soon as the system is polarized, the spin 
gap closes and all the above-mentioned correlation func- 
tions decay as power laws. Due to the unequal number 
of up and down particles, the corresponding Fermi vec- 
tors hp a = N a 7r/L = n (1 ±p) 7r/2 are different, and one 
might expect that the singlet superfluid correlations still 
remain important, albeit now composed of pairs with fi- 
nite momentum Q = — kpi\ = npn. This is the 
prediction of a bosonization analysis [27]. 

III. CORRELATION EXPONENTS FOR p > 

It is well known [28, 29] that the attractive Hubbard 
model can be mapped onto the repulsive one under a 
staggered particle-hole transformation of one of the two 
species. To be specific, let us use the transformation 

c 3"t = ( — -O^zr ' C U = ' 

The sector with particle density n and polarization p 
of the attractive model maps onto n = 1 — np and 
p = (1 — n)/(l — np) on the repulsive side (symbols with 



a bar denote quantities for the repulsive model). Un- 
der this transformation, density correlations in the at- 
tractive model map onto longitudinal spin correlations in 
the repulsive model, while singlet superfluidity is mapped 
onto transverse spin correlations, see Table I. The Fermi 
momentum of the transformed species is shifted, kp^ = 
7r — fcp,Tj while kp,i = kp,i is left unchanged. 

With this correspondence in hand, we can use the re- 
sults of Frahm and Korepin [30] on the correlation expo- 
nents for the repulsive Hubbard model in a magnetic field 
and translate the behavior of the correlation functions to 
the attractive case. In general, the leading asymptotic 
behavior of an equal-time correlator can be written as 

C (l) = A cos ( a Z k Fl l + Ptj k Fjl) (2) 

^ ' ^ |^|A € (Q!,/3) ' v ) 

where the subscript £ denotes the type of fluctuation, 
see Table I. The correlation exponents can be derived 
from the elements of the dressed charge matrix, while the 
coefficients a and (3 are constrained by the parity of the 
quantum numbers of the intermediate states [30]. Note 
that Eq. (2) only gives the leading term of the asymptotic 
behavior of the correlation functions. In principle there 
are also higher harmonics present. However, their corre- 
lation exponents being larger, one can safely neglect these 
contributions if one is interested in the asymptotic behav- 
ior. Our numerical results presented in the next section 
show that the dominant FFLO mode occurs at the ex- 
pected single momentum Q. A numerical solution of the 
Bethe-Ansatz equations for the Hubbard model yields the 
elements of the dressed charge matrix, which allow one 
to sketch the the evolution of the correlation exponents 
with density and polarization [30]. Our results for the 
dominant two-particle correlators obtained in the strong- 
coupling limit U ^> t are shown in Fig. 1(a). First, for 
p = 0, we plot the values of Ass(& = 0) for three different 
densities (filled symbols). In this limit, the CDW corre- 
lation exponent is given by Acdw(2&f) = 1/Ass(0) > 1. 
Second, for p — > + it has been predicted by Yang [27] 
that the FFLO Ass (&ft — ^fi) exponent is given by the 
SS exponent at p = plus 1/2 (open symbols). Our data 
confirm this result, see also Ref. 32. Then, upon increas- 
ing the population imbalance, this correlation exponent 
decreases monotonically until it reaches unity as p — > 1. 
The leading contribution to the asymptotic behavior of 
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FIG. 1: (Color online). Correlation exponents of the polarized attractive Hubbard model, (a) Dominant two-particle correlation 
exponents A^ (k) in the strong- coupling limit (\U\ ^> t) as a function of polarization p for different densities n. Singlet 
superfluid fluctuations (SS) at momentum |&f,t — ^f,i \ are dominant, while 2/cfj. charge density waves (CDW) are the subleading 
fluctuations. Note that the exponents Ass differ by 1/2 for p = and p — > + ; as predicted in Ref. 27. (b) Evolution of the 
Green's function correlation exponents with polarization for different densities in the strong- coupling limit. The SS exponents are 
shown again for comparison. The Green's function of the majority particles (Aq^) always decays more slowly than the minority 
one (Aqi)- (c) Phase diagram for an infinite number of weakly coupled chains as a function of density and polarization for 
various interactions U. In the FFLO phase, SS correlations are dominant, while in the Fermi liquid phase, the asymptotic 
decay of the Green's functions is slower than that of any two-particle correlator. For finite U, only the boundaries determined 
by the majority particles (Ass = Aqt/ are shown. SS denotes the conventional balanced s-wave superfluid for all n atp — O. 



the CDW correlations is found for wave vectors 2kp[, but 
their correlation exponent is always sizably larger than 
the FFLO exponent and increases towards 2 as p —> 1. 
Note that even in the limit p — » + the density correlation 
exponents at 2kp[ and 2kp^ do not become equal, sig- 
naling the absence of spin-charge separation at all p > 
in the general case n / 1 [32, 33]. 

Although longitudinal SDW and TS correlations also 
decay as power laws, they are not shown here because 
they decay more rapidly (with exponents A > 2). These 
results clearly show that FFLO pairing is the dominant 
two-particle correlation for any nonzero polarization; 
hence, there is no Chandrasekhar-Clogston limit [22, 23] 
beyond which the pairing breaks down in this strictly ID 
model. 

While the above analysis of the two-particle correlation 
exponents shows that pairing fluctuations dominate, the 
ID character of the system prevents the formation of a 
truly long-range-ordered SS phase. In higher dimensions, 
the restriction to power-law decay is lifted and one ex- 
pects that by coupling an infinite number of chains with 
a suitable interchain interaction, a long-range ordered 
FFLO phase can be stabilized. Following the standard 
analysis of the dimensional crossover, see, e.g., Ref. 31, 
we compare the dominant two-particle correlations with 
the Green's functions and argue that upon lowering the 
temperature of the system, one either enters a phase 
dominated by the physics of ID, if the two-particle corre- 
lations decay more slowly than the Green's function, or, 
one flows first to a higher dimensional Fermi liquid if the 
opposite condition is satisfied. 

In Fig. 1(b), we plot the dependence of the SS and the 
two Green's function correlation exponents on the po- 



larization p for different fillings n in the strong-coupling 
limit U ^> t. In the unpolarized system, the Green's func- 
tion decays exponentially because of the spin gap. In the 
polarized regime, the exponent of the majority particles 
is smaller than the minority one, but both tend to 3/2 
in the limit p — > 1~. Fig. 1(c) shows the phase diagram 
obtained by comparison of the single- and two-particle 
correlation exponents as a function of density and polar- 
ization for various interactions U . Here we go beyond 
Ref. 27, which neglected the fact that spin and charge 
are coupled immediately for all p > sufficiently away 
from half filling n = 1. Because of the difference between 
the majority and minority exponents, the phase bound- 
ary between the FFLO and the Fermi liquid regime is 
smeared out, i.e., there are two boundaries determined 
by Ass = and Ass = Aqj,. For simplicity, the two 
boundaries are only shown for the strong-coupling limit, 
while for finite [/, only the crossover with the majority 
particles (Ass = Aqt) is indicated. In dilute systems, 
FFLO pairing is observed for all polarizations, while at 
higher fillings, the superfluid phase is only established in 
situations with strong population imbalance. This leads 
to two disconnected patches of FFLO ordering for inter- 
actions U < 4, as can be seen in Fig. 1(c). For instance, 
for U = —1, the FFLO phase is realized in the region 
n < 0.1, but there is also a small reentrant patch in the 
vicinity of the fully polarized half-filled system at n ~ 1 
and p ~ 1. Note that this simple comparison of the cor- 
relation exponents does not take into account residual in- 
teractions that might destabilize the Fermi liquid phase 
at low temperatures and neglects the precise nature of 
the interchain couplings. For a parallel discussion of the 
dimensional crossover, we refer the reader to a very re- 
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cent work by Zhao and Liu [34] , in which a similar phase 
diagram has been presented. 

IV. NOISE CORRELATIONS 

Recently, Altman et al [24] suggested using spatial 
noise correlations as a universal probe of fluctuations in 
cold atomic gases. In experiments, these correlations can 
be extracted from a statistical analysis of many time-of- 
flight images of the expanding cloud of atoms, once they 
have been released from the trap. A key ingredient of the 
method is the fact that the spatial density distribution 
of the expanding cloud becomes proportional to the mo- 
mentum distribution in the interacting system after the 
trap is released. For the lattice model we consider in this 
work, spatial noise correlations are defined as 

G aa > (k, k') = (rik^rik',*') - (n k ,a) {n k >,<T>) , (3) 

where n^^ — c\ a c k a is the occupation number operator 
in momentum space. We will also refer to the total noise 
correlation 

G(k,k') = J2 G aa ,(k,k') . 

A. Strong-coupling limit 

We now proceed to a quantitative numerical evalu- 
ation of the noise correlations in the attractive Hub- 
bard model with strong interaction U/t = — 10, at filling 
n = 3/4 for different polarizations using a coordinate- 
space-based DMRG algorithm with open boundary con- 
ditions [35, 36]. For an in-depth discussion of the method 
and the noise correlations of various phases of extended 
Hubbard models with balanced spin populations, we re- 
fer the reader to Ref. 37. To determine the noise corre- 
lations of a system with L = 64 sites, one has to calcu- 
late 4L 4 « 67 • 10 6 four-point correlators. This number 
could be reduced by exploiting symmetries; however, we 
choose not to do so, to get unbiased results for the sum 
rules [37], which we use to check convergence. While a 
system with 64 sites seems rather small for DMRG calcu- 
lations, these sizes are, however, sufficient for the purpose 
of the present work. The main features of the noise corre- 
lations are already well established in these systems and 
finite-size scalings do not show any ambiguities. For such 
moderate system sizes, keeping up to 400 states results 
in a discarded weight smaller than 10 _T . 

Let us start the analysis with the unpolarized system, 
whose shot noise is shown in Fig. 2 [row (a)]. In this 
case, fermions tend to pair in on-site singlets, leading to 
a gap in the spin sector and thus exponentially decaying 
SDW and TS correlations. Since the correlation expo- 
nents Ass < ^cdw, the SS correlations are dominant 
over the CDW correlations. This can be seen in the shot 
noise. Quite generally, particle- hole fluctuations appear 



as negative signals along k' = k±q, where q is the order- 
ing wave vector, while particle-particle correlations give a 
positive contribution along k' — —k±q [24, 37, 38]. The 
dominant positive signal of the total shot noise G around 
(±kp, =F&f) 5 shown in Fig. 2(a,i), thus clearly indicates 
that SS correlations dominate over CDW fluctuations. 

As soon as the system is polarized, one recognizes two 
important modifications of the noise correlations with re- 
spect to the unpolarized case, see Figs. 2 [rows (b) to (f)]. 
Most noticeably, the distinct comb present in the ||- 
channel splits into two ridges peaked at (±fcpf, T^F|) 5 
which move away from the anti-diagonal with increasing 
imbalance. The splitting Q = \kp^ — kpi\ between them 
directly corresponds to the momentum of the Cooper 
pairs. This is the fingerprint of the FFLO state. The 
abrupt change from a comb with roughly constant height 
to pronounced peaks is due to the sudden increase of the 
correlation exponent upon entering the polarized phase. 
In the same spirit, we note that the peaks become more 
pronounced for increasing p, in agreement with the fact 
that the FFLO exponent decreases with increasing p, see 
Fig. 1(a). Secondly, traces of the Fermi surface, which 
are completely absent in the unpolarized system due to 
the presence of a spin gap, appear in all parts of the 
noise correlations as a nine-tile pattern best visible in 
Figs. 2 [column (ii)]. These Fermi surfaces are charac- 
teristic of a system with gapless excitations [37]. Based 
on the fact that particle-hole correlations give rise to a 
negative signal in the noise correlations while particle- 
particle contributions are positive, we expect four of the 
nine tiles to have negative values. This is in agreement 
with our numerical results. A density plot, along with 
plots of slices along the diagonal of Gjj, allows one to 
study the FFLO phase in more detail. First, as de- 
rived before, the peaks of the SS signal in the | j-channel 
are located at (ifc^ , =F&F|) and the comb extends along 
k' = —k ± (kpi — kpi)- In the finite-size scalings of the 
SS signals parallel to the ant i- diagonal, see Figs. 2 [col- 
umn (hi)], one can see that the correlations are peaked 
at (±/ci?t,T^F|) and drop down to a constant positive 
background rather quickly away from this point. This 
phenomenon reflects the fact that the pairing of particles 
is stronger the closer the two momenta are to the cor- 
responding Fermi vectors. Going along the direction of 
one of the arrows in the density plots, see Figs. 2 [column 
(ii)], one encounters a first peak representing the pairing 
of down particles at momentum k' « kp[ with up par- 
ticles at momentum k « — fcpf. Continuing along this 
line, the k' momentum crosses the horizontal axis and 
approaches the Fermi point —kpi- Although pairing in 
this region is kinematically possible, it is clear that down 
particles around k' « — kp[ are preferably matched with 
up particles from the other side of the Brillouin zone at 
kp"\- This illustrates the fact that the pairing correla- 
tions are closer to a Larkin-Ovchinnikov [4] type involv- 
ing both ±Q pairs than a Fulde-Ferrell [3] type with a 
single Q. Concerning experiments, it is encouraging to 
see that state-selective measurements of the shot noise, 
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FIG. 2: (Color online). Noise correlations in the attractive ID Hubbard model for several spin polarizations p at density 
n = 3/4. In the unpolarized case [row (a)], one identifies dominating s-wave superfluid correlations in G^i, (a,ii), along the 
anti- diagonal, indicating the presence of pairs with zero total momentum. In a population-imbalanced system [rows (b) to (f)], 
the FFLO state is clearly visible in G^i [column (ii)] and also G [column (%)] as strongly peaked signals that shift apart as the 
polarization increases. The separation Q [indicated in (f,ii)] between the peaks is equal to the total momentum of the pairs. 
Finite- size- scaling extrapolations of LG^i taken along q — —q — fcpt + &F|, indicated by the arrows shown in column (iii), 
illustrate the growth of the signal with increasing system size. The momentum distributions [column (iv)] show a clear difference 
between the gapped (p = 0) and the gapless (p > 0) phases, but do not contain direct evidence for the presence of an FFLO 
state. 
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(iii) U=-4 (iv) U = -10 




FIG. 3: (Color online). Off-diagonal noise correlations G^i 
in an unpolarized system with L = 64 sites at filling n — 
3/4 for different interactions U. The Fermi surfaces (nine- 
tile pattern), visible for small interactions, become weaker and 
finally disappear for strong couplings. SS correlations develop 
around opposite Fermi points and start to spread out along 
the whole anti-diagonal with increasing attraction. 



i.e., separate measurements of the different G ac7 > are not 
necessary in order to identify the FFLO phase, since at 
least in the highly polarized regime, the superfluid fea- 
tures present in the total shot noise G are well separated 
and almost as pronounced as the main signal along the 
diagonal, see Figs. 2 [bottom of column (i)]. 

The occupation number n a (k) for the two species sep- 
arately as well as their sum is shown in Figs. 2 [column 
(iv)]. The main characteristics for all nonzero polariza- 
tions is the sharp transition from empty to occupied or- 
bitals and vice- versa at momenta |fcp,j|, similar to the 
abrupt change observed in the momentum distribution 
of a metallic Tomonaga-Luttinger liquid, see Fig. 6(d) in 
Ref. 37. In contrast, in the spin-gapped phase at zero po- 
larization, the Fermi points are invisible. We note that 
these density profiles do not contain direct evidence for 
the presence of dominant FFLO correlations, since, for 
example, the population-imbalanced repulsive Hubbard 
chain would show a very similar n a (k) without having 
dominant FFLO correlations. 



B. Away from strong-coupling 

So far, we have focused on the strong-coupling regime, 
taking U/t = —10. Reducing the attractive interaction 
allows one to study the influence of the correlation ex- 
ponent on the shot noise. Fig. 3 shows the off-diagonal 
correlator G^ of the unpolarized system for interactions 



U/t = -1, -2, -4, and -10 at filling n = 3/4. According 
to Ref. 31, these values of the interaction approximately 
correspond to Luttinger exponents K p = 1.1, 1.2, 1.3, 
and 1.4. Note that SS correlations decay with a power- 
law exponent 1/K p . For small U, the two peaks are very 
narrow and strongly confined to the area surrounding op- 
posite Fermi points. With increasing J7, the interaction 
can scatter particles from well below the Fermi surface, 
which leads to the spreading out of the base of the peaks. 
For strong interactions, the Fermi points are barely rec- 
ognizable, the ridge extending almost uniformly along 
the whole anti-diagonal. It is thus the width of a peak 
or dip rather than its height that allows one to make a 
statement about the strength of the interaction and thus 
the correlation exponents. Noting that widths and mag- 
nitudes are related by a sum rule that ensures vanishing 
volume under the shot noise [37]; it is clear that with 
increasing width, the height of the peaks must decrease, 
as can be seen in Fig. 3. 



C. Effects of a harmonic trap 

In experiments, the atomic cloud is confined in a har- 
monic trap of the form 

L 2 
fftrap^EO"^) ( 4 ) 

1=1 ^ ^ 

with a trapping potential V. Several groups [39, 40, 41] 
have recently analyzed the FFLO pairing in ID systems 
subject to a confining potential. Compared to the ho- 
mogeneous case, the FFLO phase does not extend over 
the whole chain, but is confined to an interior part of the 
system: For small densities, an FFLO state is formed in 
the center of the system, surrounded by fully polarized 
wings of unpaired majority particles [40, 41]. In con- 
trast, for denser fillings, i.e., n ~ 1, the fully polarized 
wings remain, but the FFLO phase is pushed away from 
the center of the trap and partly replaced by a polar- 
ized metallic phase of freely moving minority particles 
in a uniform background of majority fermions [40]. In 
order to detect an FFLO signal in the noise correlations 
under realistic experimental conditions, it is thus advan- 
tageous to work in the low-density region. To study the 
effects of a harmonic potential in this regime, we have 
calculated the local densities for systems with L = 64 
sites in the strong-coupling limit U/t = —10 at filling 
n = 1/2 and polarization p = 1/2 for various depths of 
the trap, see Figs. 4[row (a)]. Even for shallow confine- 
ment [Fig. 4(ai)], the polarized wings are already visible. 
It is clear that by pushing some of the majority particles 
away from the central region of the trap, one changes the 
density and the polarization in the FFLO region. It is 
useful to introduce an effective density n e ^ = N e ^/L e ^ 
and polarization p e R = AN e ^/N e ^ at the center of the 
trap, where the density of the minority spins is greater 
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(i) V=0 .000 1 (ii) V=0 .0005 (iii) V=0 .00 1 (iv) V=0 .002 




FIG. 4: (Color online). Local densities for systems with L = 64 sites in the strong- coupling limit U/t — —10 at filling n — 1/2 
and polarization p = 1/2 /or various depths of the trap V [row (a)]. Even for small confinement, some majority particles are 
pushed away from the center and are arranged in fully polarized wings. In the central area of the trap (shaded area), an FFLO 
state at an effective filling n e jj and polarization p e jj is formed, as can be deduced from the characteristic signal in the off-diagonal 
noise correlator [row (b)] and in the total shot noise G [row (c)]. 



than some small arbitrary cutoff e = 10 3 , by defining 

i 

AiV eff ^ 0(n id - e) (n i)t - n id ) , 

i 

together with the effective length of the core 
^eff = ^20(n id - e) , 

denoting the usual step function. As illustrated in 
Figs. 4[row (a)], the effective density in the center of the 
trap increases with increasing confinement, while at the 
same time the effective polarization is reduced. The pres- 
ence of FFLO pairing correlations in the central region 
of the trap can be deduced from the characteristic signal 
present in the noise correlations G^ shown in Figs. 4 [row 
(b)]. In fact, one can convince oneself that the separa- 
tion between the two ridges corresponds very accurately 
to the momentum of the pairs in the middle of the trap 
Qeff = ^eff Peff 7T- For a shallow trap, the effect is only 
subtle, because an increase in density is accompanied by 
a reduction of the polarization. Both the increased den- 
sity as well as the reduced polarization are unfavorable 



to the visibility of the FFLO signal in the noise corre- 
lations. First, the correlations decay faster than in the 
homogeneous system, see Fig. 1(a), and, second, because 
the signal is proportional to the number of paired parti- 
cles, there is a further reduction of the FFLO signal due 
to the smaller effective volume v G s — N^jN in which 
pairing correlations dominate. Despite these adverse ef- 
fects, the total shot noise shown in Figs. 4 [row (c)] clearly 
illustrates that, for moderate traps, the FFLO signal is 
still detectable. 



V. CONCLUSION 

To conclude, we have derived the correlation exponents 
of the ID attractive Hubbard model at finite spin po- 
larizations and have shown that the FFLO type pairing 
correlations are dominant for polarizations p > 0. Fur- 
thermore, we have analyzed the spatial noise correlations 
(density-density correlations in momentum space) in 
view of experimental realizations with ultra-cold atoms. 
We have found that the ID analog of the inhomogeneous 
FFLO state leaves a distinct fingerprint in the shot noise. 
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